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The Bose-Hubbard model of a two-fold degenerate Bose gas is studied in an optical lattice with 
one particle per site and virtual tunneling to empty and doubly-occupied sites. An effective Hamil- 
tonian for this system is derived within a continued-fraction approach. The ground state of the 
effective model is studied in mean-field approximation for a modulated optical lattice. A dimerized 
mean-field state gives a Mott insulator whereas the lattice without modulations develops long-range 
correlated phase fluctuations due to a Goldstone mode. This result is discussed in comparison with 
the superfluid and the Mott-insulating state of a single-component hard-core Bose. 

I. INTRODUCTION 

An ultracold Bose gas is brought into an optical lattice, created by a stationary Laser field-'^. If the corresponding 
periodic potential is sufficiently strong, such that tunneling of atoms between the potential wells is strongly suppressed, 
the phase coherence of the Bose gas is destroyed due to the repulsive interaction between the bosons. In this case the 
Bose gas becomes a Mott insulator. This ground state is characterized by a fixed number of bosons in each lattice well 
and strong incoherent phase fluctuations of the quantum state, in contrast to the phase coherence of the Bose-Einstein 
condensate with a strongly fluctuating local particle number. The Mott insulator and the transition to the superfluid 
state were discussed theoretically some time agc^''^''* and observed experimentally recentljA. 

Besides the phase of the bosons and the local particle number there can be other degrees of freedom in an ultracold 
gas of bosonic atoms which may also establish some long-range ordering. A possible candidate for such a consideration 
is a Bose gas in an optical or a magnetic trapS^, where the fluctuations between nearly degenerate hyperflne states 
represent an additional degree of freedom. This can play a role in establishing new types of ordering, similar to the 
spin degree of freedom in fermionic systems. Here the case of a strongly interacting Bose gas with one particle per 
site will be considered. Using a Bose-Hubbard model, the interaction and the chemical potential of a grand canonical 
ensemble are ajusted such that there is one particle per optical lattice site. According to the statements given above, 
the fixed number of particles per site would represent a Mott insulator. On the other hand, the local particle number 
TT-r is a sum of the particle numbers of both components (represented by a "spin" T or j) 

fir — rir^i + nr,i- 

The individual particle numbers of the two components are fluctuating quantities and can lead to a new state when 
long-range correlations develop. A similar situation can be found in a single-component hard-core Bose gas in an 
optical lattice. Then each lattice site is either empty with the quantum state |0) or singly occupied with Formally 
these two states correspond with the states | t) and | |) of the two-component gas. It is known from analytio^''^ and 
numerical calculations^ that the hard-core Bose gas develops a superfluid phase for abitrarily small tunneling rates if 
the states |0) and |1) are degenerate, i.e. when (rir) = 1/2. This behavior will be discussed in Sect. 3.1. Using the 
formal correspondence with the two-component Bose gas the development of a long-range correlated state would also 
take place for any tunneling rate if the states | t) and | J,) are degenerate. A consequence would be that these two 
states could easily separate in space, leading to an entangled state in the Bose gas. 

Ordering phenomena can be studied perturbatively, starting from isolated potential wells of the optical lattice and 
systematically turning on the tunneling between these potential wells. Because of the degeneracy with respect to 
the spin degrees of freedom in the isolated wells this requires a degenerate perturbation theory. Instead of using a 
perturbation theory for the tunneling Hamiltonian a continued-fraction approach will be applied in the following. 

The central aim of this paper is to derive an effective Hamiltonian for the two-component Bose gas with (rir) = 1 and 
to discuss the properties of this system. In this case a particle is a superposition of the states | t) and | |). Without 
tunneling between lattice sites it represents a "paramagnetic" state, i.e. a state without ordering. Tunneling, on 
the other hand, preserves the "spin": starting with a boson of a given spin this boson will spread to neighboring 
potential wells. This process requires virtual states which are empty or occupied by more than one particle. Here only 
occupation with two particles will be allowed to keep the calculational effort low. However, the continued-fraction 
approach proposed in this paper can be extended to higher orders of occupation. 

After introducing the two-component Bose-Hubbard model in Sect. 2 the single-component hard-core Bose gas is 
considered for comparison in Sect. 3. and treated in mean-field approximation in Sect. 3.1. The ground states of a 
two-component Bose-Hubbard model with (rir) < 1 and a single-component hard-core Bose gas are evaluated for a 
two-site system in Sect. 3.2. Then the effective Hamiltonian for a projected system of interacting particles is derived 
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by a truncated continued fraction in Sect. 4 and applied to the two-component Bose gas with (rir) = 1 in Sect. 5. 
Finally, the two-site system (Sect. 5.1) and a mean-field approximation (Sect. 5.2) are discussed for the effective 
model. 

II. THE MODEL: INTERACTING BOSE GAS 

In order to describe a multi-component interacting Bose gas the interaction of bosons can either be included as a 
hard-core interaction or in terms of the Bose-Hubbard model. Crucial is only that a repulsive interaction stabilizes 
an incompletely filled lattice with vanishing compressibility. From this point of view the actual form of Hq is not 
important except for the fact that it depends only on the local particle number at site r 

with boson creation (annihilation) operators (a). For a more specific discussion, a Bose-Hubbard model with 

Ho = Y^i-IJ-rir + Ufiririr - 1)] (1) 

r 

shall be considered, where /x is the chemical potential and U > the interaction constant. Eigenvalues of the local 
particle number are n = 0, 1, ... with corresponding energies per lattice site 

i<^(n) = — /xn -I- Un{n — 1). 

For < II < U , the case considered throughout this paper, the lowest energy is E(l) = —jjl for n = 1, and next higher 
energies are E{0) — and E{2) — 2{U — /i) and even higher energies for n > 2. 
The dynamics of the bosons is described by the tunneling Hamiltonian 

Hl = - ^ Tr,r' ^ taal^„ar',a, (2) 
<r,r'> <T=l,1 

where < r, r' > arc nearest-neighbor sites on the optical lattice. The parameter < Tr.r' < 1 describes a modulation 
of the optical lattice. The interaction is crucial for finding new physical states, since the noninteracting Bose gas (i.e. 
for [7 = 0) the Hamiltonian Hi gives only two independent Bose-Einstein condensates. 

A grand-canonical ensemble of bosons at the inverse temperature defined by the partition function 

Z = Trexp{-f3H), (3) 

can be used to evaluate the average density of particles as 

1 dlogZ 

71 = 

(3N dn 

with the number of lattice site N. Without tunneling (i.e. for = 0) the average density of particles n gives 
dn/dfj, = (i.e. incompressible states) for all non-integer values of n/U. 

III. SINGLE-COMPONENT HARD-CORE BOSE GAS: SUPERFLUID AND MOTT-INSULATING 

STATES 

Before discussing the two-component Bose gas the single-component hard-core Bose gas shall be considered because 
its properties are known from a number of other approaches. The hard-core interaction can be described by a bosonic 
creation (annihilation) operator (A) with the additional condition = 0. The Hamiltonian of the hard-core 
Bose gas then reads 



Hhcb = -t ^ AlAr' - n'^AlAr 

<r,r'> r 



(4) 
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and acts on the Hilbert space with |0) and |1) as basis states at each lattice site. There is a close formal connection 
between hard-core Bose and spin-1/2 operators, since one can write 

S'-" ^ {A + A'^)/2, Sy ^i{A- A^)/2, S'- = A- 1/2. (5) 

The Hamiltonian can be expressed in terms of these spin operators as an XY model with a magnetic field in z direction: 

<r,r'> r 

For sufficiently large values of the ground state is ferromagnetic: if /i > the magnetization is {S^) > and vice 
versa. Therefore, for /i > (/i < 0) these ground states correspond in terms of the hard-core bosons with a completely 
filled (empty) lattice, where the filled lattice represents a Mott-insulating statc^'"^'^. /i = is a marginal situation, 
where the bosons develop a superfluid state for any positive value of the tunneling rate t. More general, the superfluid 
state persists if the tunneling dominates. 

The advantage of the spin representation is that it provides a simple qualitative picture for the existence of a 
Mott-insulating state and a transition to a superfluid state. A more quantitive description is obtain from a mean-fleld 
approximation which is discussed in the next section. 



A. Mean-field approximation of the hard-core Bose gas 

A possible complex mean-field state for the hard-core Bose gas is 

l*A./F> = n [e'"- cos(?7r) + e^'^- sin(r,r)4] |0) (6) 

r 

from which matrix elements can be calculated. For instance, the tunneling term in the Hamiltonian Hhcb gives 

(^MfI^JA' + Aj,Ar|^'MF) = 2c0s(ar - Or' ) COs(77r ) sin(77r ) C0s(77r' ) sin(77r' ) , 

where the phases appear only in the phase difference — (pr ~ "fpr- Thus this term of the Hamiltonian has a global 
U{1) symmetry because it is invariant under a shift ip^ ipr + and ■(/'r — > + The other matrix element of 
the Hamiltonian Hhcb is independent of the phases: 

(*MF|^];^r|*J\/F> =sin2?7r. 

The mean-field expectation of the Hamiltonian with the homogeneous mean field rj reads 
{'^MF\Hhcb\'i'MF) = - sin^ ?7 t{l - sin^ ry) ^ cos(ar - ar') + X] ' 

<r,r'> r 

The ground state is cos(ar — ar') = 1 and 

r for ^ < -2dt 

sin^ r]= I 1/2 + fi/Mt for -2dt < fi< 2dt . (7) 
I 1 for 2dt < /i 

The expression 

(*MF|Ar|*MF> = e*('^^-'^-) cosr,sin77 (8) 

is an order parameter for a superfluid state which vanishes in the cases sin^ rj = 0,1. It should be noticed that in this 
regime (5'a/f|^^/ic&|^j\/f) is independent of Or- This reflects the absence of a superfluid. In the case sin^ = it is 
an empty state and for sin^ 77 = 1 a Mott insulator with one particle per site. This result is in agreement with Monte 
Carlo simulations^. 
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B. Hard-core vs. Bose-Hubbard model: a two-site system 



To study the tunneUng between neighboring potential wells of the optical lattice a two site model with = = t 
is considered. This was already discussed in great detail for the model under consideration, using the leading order 
of an expansion in t'^/UL^. Here the study shall be performed without assuming t'^/U <C 1. As a first example the 
single-component hard-core Bose gas is considered in terms of the Hamiltonian in Eq. Q . The ground states are 

nO,0) for ^ < 

l^-o) = < (|0, 1) + |1, 0))/^/2 for -t < /i < i • 

for^>t 

This result corresponds with the three different mean- field ground states in 0. Thus already the two-site system 
indicates the three phases of the full d-dimensional lattice: the empty lattice, the condensate and the n = 1 Mott 
insulator. An analogous calculation for the two-component Bose-Hubbard model with empty and singly-occupied 
sites gives the ground states 

( |0, 0) for n<-t 

l*o) = h|0,t) + | T,0))/V2, (|0,i) + U,0))/V2 ior-t<fi<t . (9) 

Ui,i>, U,T>, IT,i>, IT,T> for/i>i 

For —t<i.j,<t there are two (degenerate) ferromagnetic states, where the degeneracy can be lifted by an infinitesimal 
magnetic field. It is expected that these are the states which form a condensate. The four-fold degeneracy for /i > t 
may be lifted when a virtual tunneling through empty and doubly-occupied sites is included. This degeneracy also 
raises the question whether or not an analogue of the superfluid ground state is allowed due to virtual tunneling 
in this regime. A reason for having long-range correlations is that one of the two components, e.g. | |), can be 
formally considered as an empty site, the other component as a hard-core boson and since the hard-core Bose gas has 
a superfluid state for sufficiently large tunneling rate. To study this regime, a projection of the trace in the partition 
function to singly-occupied states is considered subsequently. This projection allows virtual tunneling through empty 
and doubly-occupied sites. 



IV. A CONTINUED-FRACTION APPROACH TO THE PROJECTED PARTITION FUNCTION 

The many-body system is defined by the Hamiltonian H and the transfer matrix . Physical quantities at inverse 
temperature /3 are derived from the partition function Z defined in Eq. (Q. A continued-fraction approach shall be 
developed in this section to derive an effective Hamiltonian from H that describes the physics of a projected transfer 
matrix Pqc^^^ Pq. Although H is bounded from below it may have negative eigenvalues. A positive operator can be 
obtained by adding a constant diagonal term E to shift the ground state energy Eq to positive values £'(1) = Eo + E. 
Then the transfer matrix can be represented by the integral 



Thermodynamic properties at low temperatures are dominated by the ground states and low-energy excitations. 
The trace of the grand-canonical partition function Z includes states with all possible number of bosons. For the 
Hamiltonian iJ the highest statistical weight comes from the ground state, li H = Hq + Hi with a perturbation Hi, 
the Hilbert space is projected on the degenerate ground states of the Hamiltonian Hf) by Pq. In the example of Sect. 
2 the Hamiltonian preserves the number of particles. Therefore, in this case it is expected that the system with fixed 
particle filling gives the dominant contribution to the trace, especially at low temperatures. This situation can be 
described by the Po-projected partition function with the corresponding trace Ttq 

PE roc 

Z = Tro{Poe-^"Po) = — / e'^'Tro[Po{z - i{H + E))-^Po]dz. (10) 

Assuming that H is implicitely shifted by E, the Po-projection of the resolvent [z — iH)~^ reads 

Po(z - iH)-^Po = (Poiz - iH)Po + PoHPiiz - iH)^^ PiH P^)^^ (11) 

with Pi = 1 — Pq. \ is the inverse with respect to the Po,i -projected space. This identity can be directly shown by 
a multiplication of the matrix and its inverse. It can be generalized to a recurrence relation (s. Eq. I|23|) of Appendix 
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A) if the Hamiltonian H satisfies the special conditions (|22|) . For the two-component Bose gas this is indeed the case, 
since the Hamiltonian H — Hq + Hi of Sect. 2 is of the special form 

TT _ ( PoHqPo PoHiPi\ . ^ 



PiHiPo PiHPi 

where Pq is the projection to the degenerate ground state of Hq with Q < fi < U, i.e. one particle per site. Thus the 
Hamiltonian Hi is responsible for an interaction between the Pq- and the Pi-projected Hilbert spaces. Hq, on the 
other hand, acts only inside the Po-projected space. Starting from singly-occupied states, the tunneling Hamiltonian 
Hi can only create a pair of an empty and a doubly-occupied site (PEDS). If P2 is the projection from Pi to one with 
only a single PEDS, the second term in the inverse matrix of Eq. Ullfl reads 

PoHiPiiz - iH)i'PiHiPa = PoHiP2{z - iH)i'P2HiPo. (13) 

In general, the operator P2k+iHiP2k creates a new PEDS in the Hilbert space with k PEDSs. Thus the continued- 
fraction representation of Appendix A, applied to the two-component Bose gas at single-particle filling, is based on 
the creation of PEDSs. 

In order to truncate the continued fraction the creation of new PEDS and multiply-occupied sites is excluded. This 
is related to the approximation 

P2iz - ^H)l^P2 « P.(z - ^H,)l^P2 ^ , _ ,[(^ _ ^M^) ^ ^(q) + i^(2)] 

where E{Q) and E{2) are the energies Hq+ E for the empty and doubly occupied sites of Sect. 2. With Eqs. Ullf) 
and H13|l this gives for the Po-projected resolvent of the partition function 

Po(z - iH)-^Po « (Po(z - iH^)Po + PoHiP2[z - iHo]-^P2HiPo)^^ 



= {z- iEo + PoH^Po/[z - i[{N ~ 2)Eil) + E{0) + E{2)]])o\ 

Here it has been used that 

P0H1P2H1P0 - PoHfPo 

which follows from the fact that Hi is off-diagonal with respect to the Pq- and P2-projected Hilbert spaces. Then the 
Po-projected partition function reads for low temperatures (i.e. f3 ^ 00) (s. Appendix B) 



(1 - AEH-/f) (14) 



with the effective Hamiltonian 

Heff = -((AS)2 + Poi/2pg)i/2 = _^E{1 + PoHfPo/AE)^^^ 

and 

AE= [E{Q) + E{2)]/2^E{1). 

Since PoH^Pq is a non-negative operator, this result imphes that the ground state of —PoH^Pq is the ground state 
oiHeff. 

If PqHiPq is small in comparison with (Ai5)^, a perturbation theory with respect to the PqH^Pq can be applied 
to Hf.ff- This leads to an expansion with respect to t'^/U. In leading order the approximation is 

Heff ~ -AE - PoH^Po/i2AE), 

in agreement with the results of Refs7'^'^. 

This method of deriving an effective Hamiltonian by projecting the partition function is quite general as long as the 
Hamiltonian H has the structure shown in Eq. (|12|l . The specific case of a two-component Bose gas will be discussed 
subsequently. 



6 



V. THE EFFECTIVE HAMILTONIAN OF THE PROJECTED TWO-COMPONENT BOSE GAS 

Numerous possibilities were discussed in the literature for the creation of a two-component system in atomic 
gg^gggi0jiiji2ji3 Pqj, instance, if ®^Rb is coupled to a radiation field there are pairs of nearly degenerate hyperfine 
|F,mp.) states, namely | T> = and | i) = |2,l)i044 or | t) = |1,-1) and | i) = \2,-2)^. \ [) and | T) are 

formal notations to specify the two (almost) degenerate states. The interaction of the atoms does not depend on the 
states but only on the local density of the bosons. Therefore, the interaction in iJp of Eq. is a good description. 

In a dilute regime (i.e., for a filling less than one particle per site) the interaction is weak. This opens the opportunity 
to apply a classical approach for the two-component condensate order parameter, leading to a lattice version of the 
Gross-Pitaevskii equation. The optical lattice means that the kinetic term has a special band dispersion e(k) in 
Fourier space, depending on the lattice, instead of the fc^ dispersion in the case without an optical lattice. It is 
believed that this model has a ferromagnetic ground state with respect to the two-component bosonsi^, a result that 
is also supported by the result of the two-site system in Eq. @. On the other hand, it is well known from the theory 
of the fermionic Hubbard model that the type of spin order depends crucially on the filling of the latticeiS*ii, and can 
change, for instance, from ferromagnetic to antiferromagnetic order by changing the filling. In particular, the fermionic 
Hubbard model has an antiferromagnetic ground state at half-filling. To study this effect in the two-component Bose 
gas in an optical lattice, the strongly interacting case with single-particle filling of the lattice shall be considered here. 
Then the two-component degeneracy has to be taken fully into account and the projection approach of the previous 
Sect, for the tunneling term Hi should be applied. In this case the Hamiltonian —PqHiPq reads 

~2 ^ ■^r^r' X! t<yta'Po{al,aar',a+al,^^ar^cr)ial,j,ar'^cr'+al,^^,ar^^')Po (15) 
<r,r'> a,cr'=l,1 

and AE — U. On the Po-projected Hilbert space (i.e. the space with exactly one particle per site) the operators 

Al = Poal.^ar,iPo, A^ = Poflr^ar^T-fo 

are creation and annihilation operators of hard-core bosons, when | j) is formally identified with a vacuum state |0) 
and I t) with a one-particle state As shown in Appendix C, the operator —PoIlfPo reads in terms of the hard-core 
Bose operators as 

-PoHfPo^- J2 rl,,[^t^AlA,, + ^-l^il-AlAr)Al,A,,]. (16) 

<r,r'> 

This Hamiltonian describes a hard-core Bose gas with a repulsive nearest-neighbor interaction. The competition 
between tunneling (favors particles) and nearest-neighbor repulsion (favors a ground state with checkerboard order) 
leads to a complex situation. Similar to the single-component hard-core Bose gas, this can also be discussed in terms 
of spin-1/2 states. The representation of Hamiltonian as a spin Hamiltonian via Eq. gives an anisotropic 
Heisenberg Hamiltonian. 



A. Projected two-site model 



The situation of the projected partition function of the two-site model can be discussed and compared with the 
previous results for the hard-core Bose gas and the two-component Bose- Hubbard gas. Using the hard-core Bose 
Hamiltonian of the projected model in Eq. H16() with t-^ — = t and r^y = 1, the corresponding 4x4 matrix for 
the four differents states with | and "f at the two sites has the eigenvalues {—2t^, 0, 0, 0}. The unique ground state of 
—PqHIPq with energy —2t^ is 

l*o) = ^(U,T) + IT,i)). (17) 

The Po-projection has apparently selected a non-degenerate ground state from the four degenerate ground states of 
the Bose- Hubbard model with one particle per site in Eq. ij^. This is a consequence of the virtual tunneling to 
empty and doubly-occupied states in the model with the projected partition function, which was not included in the 
derivation of the state l^io) of Eq. I^. This state is not an eigenstate to 5*^ but has a vanishing expectation for S^. 
This reflects the fact that the ground state has no tendency to develop a ferromagnetic order. 
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The projected partition function (|14|l reads for this two-site model 



2 ^ ^-l3(2Eo + U) ^ 
J = l 

Z can be used to evaluate the average tunneling energy from 

Thus the interaction reduces the tunnehng rate of the two-site model by a factor (1 -I- U"^ /2t^)^^/'^ . 

B. Mean- field approximation of the two-component Bose gas 

Using the hard-core Bose representation of the Hamiltonian in Eq. its mean-field approximation is studied 

with the complex mean- field state (jS)). The repulsive nearest-neighbor interaction is 

(*Afi^|(l - AlA,:)Al,Ay,\^MF) = cos^r^rsin^ryr'. 

Together with the hard-core Bose Hamiltonian of Sect. 3.1 the expression in the Hamiltonian of the two-component 
Bose gas at single-filling reads in mean-field approximation 



' MFl 



{^MF\PaHlPa\^l 

2 I +2 



E2 ( / \ • • T12-2 
Tj. J./ \ t^ti C0S(ar — Cty.1 ) cos T^r Sm ?7r COS 77r' Sm 77r' H — - cos T^r Sm 77r 

<r,r'> 



The first term favors a homogeneous solution for ?7r, the second term an inhomogeneous solution, e.g., for neighboring 
sites r, r' with 

sin^77r = l, cos^?7r' = l. (18) 
In this case the first term vanishes and the remaining Hamiltonian is 



<r,r'> 

,2 , 



On the other hand, a homogeneous mean-field solution for the ground state is sin r/ = 1/2 such that 

1 ^2 _j_ j-2 

-{'i'MF\PoHfPo\^MF)h^-^ J2 Tl^'{htlCOsiar-ar') + -^-Y^y (19) 



This Hamiltonian agrees formally with the mean-field Hamiltonian of hard-core bosons in Sect. 3.1 at the point of 
degeneracy /i = 0. However, its interpretation in terms of the physical bosons, given by the Bose operators and a, 
is different. This is clearly indicated by the fact that the order parameter of a superfluid state vanishes: 

('i'MF\ar\'^MF)h = 0. 

Thus the long-range correlated phase fluctuations are not related to a superfluid state but to a spontaneously broken 
symmetry, associated with the order parameter 

{^MF\Ar\-^MF)h = (*A/i=^|a|;^^ar,T|^' 

These phase fluctuations prevent the system to become a genuine Mott insulator, since the latter is characterized by 
a gap and short-range correlated fluctuations (cf. with the Mott-insulating state of the single-component hard-core 
gas in Sect. 3.1). However, a Mott insulator can be obtained in the limit t|t| = 0. This is a limit similar to the 
Falicov-Kimball limit of the fermionic Hubbard modelii. 
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Another mean-field approximation can be constructed for a generalization of the two-site model to a modulated 
lattice model, using dimcrs with tunneling rates T^^ri = Tq as building blocks of the lattice. These dimers are weakly 
coupled with tunneling rate Tr^r' = n ^ tq- A corresponding complex mean-field state is 

1*^)= n 4(e^'^-4 + e''^-'4)|0), (20) 

<r,r'>e_D V ^ 

where D is a set of dimers {< r, r' >} with Tr ,.' — tq. This state is a lattice generalization of the two-site state of 
Eq. (|17|l . The Hamiltonian of the two-component Bose gas at single-filling (|16|l reads in this mean- field approximation 

-{^o\P^HlPo\^D)^-^^Mh E COs(^r-^r')-^^ E ^'r' ' (21) 

<r,r'>e-D <r,r'> 

For t^t^ = this result agrees with the Hamiltonian of solution H18|l but has a lower energy for any > 0. Moreover, 
the state has always a lower energy than the homogeneous mean- field state \'iMF)h- This can be summarized 

by comparing the ground-state energies: The difference between the ground-state energy of the homogeneous (Eh) 
state |^'mf)?i and the inhomogeneous (Ei) state |^j\/F)i is 

E,-E^^^^l^ ^ rl, 

<r,r'> 

and between the inhomogeneous state |^'j\/F)j and the dimerized state I^'d) is 

E^- Ed ^ -y^ E ■^r.r' • 

<r,r'>e_D 

It should be noticed that the homogeneous state \'^MF)h and the dimerized state I^'d) have the same ground-state 
energy in the case of a single-component hard-core Bose gas, provided that the dimers fill half of the lattice. This 
means that the difference of the ground-state energies in the Hamiltonian of Eq. H16|l is due to the repulsive interaction. 

The Hamiltonian of Eq. H21(l creates short-range correlated fluctuations, since the set D contains only isolated 
dimers. Therefore, it represents a Mott insulator. If the modulation of the lattice is weak (i.e. tq w ri), a summation 
of the state in Eq. H20() over different dimer configurations is required. This may lead again to long-range correlated 
phase fluctuations, since the global U{1) symmetry of the phase fluctuations can be spontaneously broken. Thus a 
phase transition from a Mott insulator at strong modulation to a state with long-range correlated phase fluctuations 
at weak modulation is expected. 



VI. SUMMARY 

A two-component Bose gas with creation operators a|, aj in an optical lattice with lattice modulations and one 
particle per site is studied. By allowing only virtual tunneling to empty and doubly-occupied sites, an effective 
Hamiltonian is derived for hard-core bosons, defined by the creation operator 

where Pq is the projector on one-particle states. The effective Hamiltonian decribes tunneling and a repulsive nearest- 
neighbor interaction between hard-core bosons. It is studied in terms of two types of mean-field states: a product of 
single-particle states and a dimerized state. The ground-state energies of a homogeneous single-particle product state 
(Eh), of a inhomogeneous single-particle product state (Ei), and of a dimerized state (Ed) are related as 

Ed < Ei < Efi, 

where the first equality sign holds for t^t^ =0 and the second for = ^i- 

The repulsive nearest-neighbor interaction prefers the dimerized state, indicating a Mott insulator at least in the 
presence of a lattice modulation. This state is characterized by short-range correlated phase fluctuations. For small 
or even vanishing modulation, however, a superposition of different dimerized states may lead to a spontaneously 
broken U{1) symmetry of the phase fluctuations. This would be accompanied by a Goldstone mode with long-range 
correlated phase fluctuations, indicating the destruction of the Mott insulator and the creation of an ordered state in 
terms of the two components of the Bose gas. 
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A similar model with N components and hard-core interaction was studied in the N ^ oo limit It has a 
Mott-insulating phase with h — I and indicates a symmetry-breaking phase for < n < 1. 

This work was supported by the Deutsche Forschungsgemeinschaft through Sonderforschungsbereich 484. 
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Appendix A 

Given is a sequence of projectors Pj {j > 0), defined by 

P2k+1 - P2k-i - P2k {k > 0) 
with initial condition P_i = 1 and with the Hamiltonian H through the properties 

P2kHP2k+l = P2kHP2k+2, P2k+lHP2k = P2k+2HP2k- (22) 

With these projectors the identity Eq. (|11|) can be iterated. The first step is to replace PqHPi{z — iH)j^PiHPo with 
the right-hand side of the identity 

PoHPiiz - iH)^'P,HPo = PaHP2{z ~ iH)^'P2HPo 

such that Eq. Hll|l reads 

Poiz - iH)-^Po = (Poiz - iH)Pn + PnHP2{z - iH)^^P2HPo)o\ 
Now the expression P2{z — iH)^^P2 on the right-hand side can be rewritten by applying again Eq. Hll|l as 

P2{Z ~ iH)l^P2 = {P2{Z - lH)P2 + P2HPi{z - lH)^^P^HP2)2^ 

with P-i — Pi — P2- Moreover, application of (|22|l to the right-hand side yields 

P2{Z - iH)^^P2 = {P2{Z - lH)P2 + P2HP^{Z - iH)^^PiHP2)2^. 

Iteration of this procedure leads to the recurrence relation 

P2k{z - iH)-^_^P2k = [P2k[z - lH)P2k + P2kHP2k+2{z - lH)-^\^P2k+2H P2k)2k ■ (23) 

Together with Eq. (|nj this gives a continued-fraction representation of Po{z — iH) ^Pq- 
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Appendix B 



To evaluate the integral (|10|1 it is convenient to define 

E, = ^AE + E{1), AE = [E{0) + E{2)]/2 - E{1). 
Moreover, the spectral representation of PqHiPq is used with the eigenvalues Xj. This leads to 



with the integral 



The poles of the integrand are 



such that the integral itself is 



^' 2m J_^{z-iNE{l)){z-iNEi) + \ 



dz. 



z± 



/(A) 



N{E{1) + El) ± ^N'^{Ei - E{l)f+A\ 
e^f'^+{z+ - iNEi) ~ e'P^-{z^ - iNEi) 



The second term of the numerator dominates at large values (3: 

e'l^^-{z^ - iNEi) 



/(A) 



z+ — Z- 



Since 



z_ = t^N^iEi- Eil))^ + AX, 



(24) 



z_ - iNEi 



and 



the expression in Eq. (|24l) reads 



N{E{1) ~ El) - y/N^{Ei- Eo)^ +AX 



El - E(l) = —AE, 
^ ' N ' 



AE 



y/{AE)^ + X 



Appendix C 

It is convenient to split the summation ^, in Eq. (|15|l into a diagonal part and an off-diagonal part: 
E Tr,r'[ X! ^lPoal^a"'r,<yar',aal,,^Po + 2t^tiPoal^ar,ial, T-ar'.iPo 

<r,r'> (T=i,T 
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The projection Pq acts individually at each lattice site r, i.e., for r' ^ r one can write 

cfPo — PoO-l^cr'^r.jjPoPoO.r' ,a0.l' rrPo 

and 

With this it is possible to define operators on the Po-projected Hilbert space (i.e. the space with one particle per site) 
as 

Al = Poal-^ar.iPQ, Ar = Poal^iar,TPo- 

When I J,) is formally identified with a vacuum state and | t) with a particle, A^ {A) is a creation (annihilation) 
operator for a hard-core boson on the Po-projected Hilbert space. Moreover, it is 

Poar,aal^Po = Po- PoaLflr.^Po, (25) 

and the operators satisfy the identities 

AlAr = Poal.^ar,iPoal^^ar,^Po = Poal^^ar,^Po = Poar,ial^^Po (26) 



ArAl = Po4_j^ar,tPoaJ iflr^xPo = PoflJ j^ar,xPo = PoflrjaJ -fPo. 



(27) 



Thus AlAr is the particle number operator for the hard-core bosons. With H25(l and H26|l the Hamiltonian —PqH^Pq 
can be written in terms of the hard-core Bose operators as 



-PoH^Pa=- ^r,r'k^i4^ 



^ll^il-AlAMlA' 



